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From calculations of the high temperature series for the free energy of the two-dimensional t-J 
model we constuct series for ratios of the free energy per hole. The ratios can be extrapolated very 
accurately to low temperatures and used to investigate phase separation. Our results confirm that 
phase separation occurs only for J/t > 1.2. Also, the phase transition into the phase separated state 
has T c w 0.25J for large J/t. 



The Hubbard and t-J models, though widely used 
to investigate high temperature superconductors, remain 
controversial when doped away from one electron per site. 
The possibility that doped holes do not form a uniform 
phase but instead phase separate into distinct high and 
low density regions on the lattice is an important issue 
that has proved difficult to settle Jl|-0. Phase separa- 
tion for physical choices of model parameters would imply 
more complicated models of 2D strongly correlated elec- 
trons are required to describe high temperature super- 
conductors. Stability of a uniform density phase would 
leave open the possibility that simple models contain the 
relevant physics without additional terms. 

While experiments have clearly observed phase sepa- 
ration in a few high-T c systems, notably oxygen over- 
doped La2Cu04+a with mobile interstitial oxygen atoms 
H , phase separation does not seem to be a universal fea- 
ture of the cuprates. However, the mechanism of phase 
separation causes holes to feel a net attraction, a possible 
precursor for the formation of stripe phases or supercon- 
ductivity. Finding an attractive interaction for holes in 
models that have predominantly strong repulsive inter- 
actions is not easy, and all known possibilities deserve 
thorough investigation. 

To investigate the properties of phase separation we 
have calculated the high temperature series for the 2D t- 
J model free energy to 10th order in inverse temperature. 
The Hamiltonian for the t-J model is 



H = -tJ2 
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where the sums are over pairs of nearest neighbor sites 
and the Hilbert space is restricted to states with no dou- 
bly occupied sites. The series is generated for a 2D square 
lattice. 

To determine the stability of the uniform phase we 
would like to investigate the ground state energy per hole 
given by 

E (S)-E£F 



e(S) 
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introduced by Emery, Kivelson and Lin [Q. Here E (S) 
is the ground state energy per site of the uniform phase 



for hole doping <5 and Eq = —1.1 6944 J is the ground 
state energy per site for the Heisenberg model H|| where 
(5 = 0. If e{8) is a monotonically increasing function 
of 5 the uniform phase is stable. If e(S) is constant or 
decreasing for a range of dopings the uniform phase is 
unstable for those values of 8. 

There are two main difficulties encountered in calcu- 
lating e(<5) from numerical measurements (exact diago- 
nalization, quantum Monte Carlo or Green's function 
Monte Carlo) of Eq(S). The first is that e{5) requires 
the subtraction of two large numbers Eq(S) and Eq F to 
determine a small number which is then divided by 5, an- 
other small number. Given statistical uncertainty in nu- 
merically determining Eo(S) (Eq F is essentially exact in 
comparison) this is a difficult task, especially for S -C 1. 
The second difficulty is that numerical calculations are 
done on small clusters. Systematic errors in Eq(S) are 
tough to estimate without knowing the finite size scaling 
of the data and whether the cluster sizes considered are 
large enough to be in the scaling limit g}|7j]. In addition 
to these difficulties, phase separation is favored on small 
clusters for S <C 1. The reduction in ground state energy 
due to the kinetic energy of the holes, which disfavors 
phase separation, is not as large on a small cluster as it is 
for an infinite lattice. On a small cluster the electron sys- 
tem reduces its energy more through local interactions, 
which for the t-J model are attractive interactions for 
antiparallel spins due to the J term in the Hamiltonian. 

High temperature series provide a means to avoid these 
difficulties. We generalize e(8) to T > by 



f(S,T) = 



F(S,T)-F AF (T) 



(3) 



where we have replaced the ground state energy per site 
by the free energy per site and limx^o f(5,T) = e{5). 
This replaces the difficulties mentioned above by the need 
to analytically continue the series to low temperatures. 
For J/t & 1 and S <C 1 we find ratios f{b~2)/ f{b~i) for two 
closely spaced dopings 5i and 62 are the best quantities 
to extrapolate. Series for ratios can be calculated exactly 
from the series for F, avoiding the need to subtract two 
large approximate numbers. The series coefficients are 
also exact for an infinite lattice so we have no explicit 
finite size effects. The ratios are extrapolated using stan- 



1 



dard Pade approximants, but only after the exact series 
for a given ratio is calculated. The doping spacing we 
use is AS = 0.025. 

By extrapolating /(^V/^i) to T — we obtain es- 
timates for e{52) / e(5\) in the uniform phase. Since high 
temperature series start at infinite temperature and only 
have information for the phase above a nonzero T c , all 
of our results are for the uniform phase. A description 
of what happens if we try to extrapolate below T c > 
is given below. Results for a range of dopings and J/t 
values are shown in Fig. 1. For the parameters consid- 
ered here e(5) < so that if <5 2 > Si and the system 
phase separates we should find e(5 2 )/e(<5i) > 1 if T c > 
or e(#2)/e(<5i) = 1 if T c = 0. If the uniform phase is 
stable we have efa) / e(S\) < 1. The 2D t-J model phase 
separates into a phase with (5 = and a doped phase 
with 5 = Sq. For phase separation we therefore expect 
e(<5)/e(0.01) > 1 immediately upon doping. 




8 

FIG. 1. Doping dependence of the energy per hole normal- 
ized to the energy per hole at 5 — 0.01, plotted as a function 
of 5 for different values of J/t. The J/t values from top to 
bottom are J/t = 0.0, 0.2, 0.4, 0.6, 0.8, 1.0, 1.1 and 1.2. The 
monotonic decrease of the ratio with increasing S indicates 
the uniform phase is stable for J/t 1.2. 

In Fig. 1, e((5)/e(0.01) < 1 and falls monotonically 
with increasing S for all J/t shown, indicating no insta- 
bility towards phase separation in the 2D t-J model for 
J/t £ 1.2. 

In Ref. a variational argument is used to support 
the presence of phase separation for J/t < 1, A vari- 
ational phase separated state was constructed from two 
pieces occupying different parts of the lattice: a Heisen- 
berg antiferromagnet for the 8 — phase and a gas of 
spinlcss fermions for the 8 — 8q phase. The energy of 
this state is then minimized with respect to 5, giving 
E (S) = E$ F -US(1- y/BnJ/t) for the phase sepa- 
rated state and S = s/BJ/nt, where B = 1.16944/2 = 
0.58472. This energy was then compared to ground state 
energy estimates for the uniform phase found by consid- 
ering a single hole in an antiferromagnet. The energy 
for the phase separated state was found to lie below the 
unform state energy for small enough J/t, and since the 



variational energy lies above the true ground state energy 
the conclusion of Ref. Q] was that the phase separated 
state is stable. Extrapolating the result for a single hole 
to a finite density of holes assumes the energy bands re- 
main rigid, a feature not obvious for a strongly correlated 
system. 



0.00 
-0.05 
E/t -0.10 







7/r=0.01 


















I 



8 

FIG. 2. Comparison of ground state energy estimates at 
J/t — 0.01. The solid line is a variational phase separated 
state which extends to So = 0.04314. The data points are 
ground state energy estimates for a uniform state calculated 
by extrapolating the high temperature series for the free en- 
ergy. At J/t — 0.01 the variational estimate lies above the 
uniform state estimate and is thus not sufficient to show phase 
separation. 

In Fig. 2 we compare our estimates for the uniform 
ground state energy to the phase separated variational 
ground state energy at J/t = 0.01. We find that our en- 
ergies lie below the variational energy for S < Sq. Note 
that from this result we cannot conclude that the uni- 
form state is stable, but only that the variational state 
discussed in Ref. jlj] is not sufficient to show phase sepa- 
ration at J/t <C 1. 
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FIG. 3. Temperature dependence of /(0.02)//(0.01) for a 
range of J/t values. From top to bottom J/t = 0.2, 0.4, 0.6, 
0.8 and 1.0. The overall temperature dependence is fairly 
small, making the temperature extrapolations more reliable. 

In Fig. 3 we show the temperature dependence of 
/(0.02)//(0.01) for a range of J/t values. Estimating the 
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low T behavior of this function is our only approximation. 
The weak temperature dependence for the ratio leads us 
to believe our results are reliable. The general trends of 
the data shown in Figs. 1 and 3 are due to the minimum 
in Eq(S) moving to smaller 8 as J/t is increased, causing 
e(6) to decrease in magnitude faster than e(0.01), though 
for the parameters shown e{8) and e(0.01) remain nega- 
tive. For values of J/t larger than shown in Figs. 1 and 3 
the ratio /(#2)//(#i) develops a spurious pole due to the 
crossing of F AF and F(Si) at T > 0. This pole greatly 
degrades the accuracy of extrapolations of the ratios at 
lower temperatures. To investigate larger J/t we need 
another method. 

The chemical potential fi = —dF/dS provides another 
means to investigate phase separation. We typically find 
\x is more difficult to extrapolate than f(52)/f($i), with 
the error in the extrapolations for fi considerably larger 
than for the ratio. For J/t > 1.2 we do see fi(6) becoming 
quite flat for 8 <C 1, as expected for a first order phase 
transition into a phase separated state. As the tempera- 
ture is lowered, \i near the critical point (critical doping 
8 C and temperature T c ) becomes flat, giving a diverg- 
ing compressibility n at the critical point. Results for \x 
are shown in Fig. 4. The flat region found in /j,(8) can 
be used to estimate the boundary for phase separation. 
However, for larger <5 distinguishing where the flat region 
ends is difficult, leading to errors in the position of the 
phase separation boundary 
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FIG. 4. Chemical potential at T = 0.2t as a function 
of doping for a range of J/t values. From top to bottom 
J/t = 0.4, 1.0, 1.2, 1.3, 1.4, 1.6, 1.8 and 2.0. For each 
J/t > 1.3 a range of dopings exists where the chemical po- 
tential is approximately constant, as would be expected near 
a critical point. 

Further evidence of phase separation at large J/t can 
be found by directly extrapolating F(6, T) to estimate 
E (S). Fig. 5 shows results for J/t = 2.0. The char- 
acteristic signature of phase separation is the reversed 
curvature observed from 8 = to 8 w 0.45. The re- 
versed curvature of Eq (8) (giving an unphysical negative 
compressibility) results from extrapolating the high tem- 
perature uniform phase F(8, T) through the T c > phase 



transition for phase separation. If T c = we would find 
instead that E (8) became linear in 8 in the phase sepa- 
rated region. 

The reversed curvature shown in Fig. 5 indicates T c > 
0, but T c is probably quite low. An indirect estimate of 
T c can be made at large J/t, above J/t = 3.4367 where 
the 2D t-J model phase separates at all densities G(J 
into regions with 6 = and 5 = 1. Here we know Eq(8) 
for all 5, since Eq(8) is the linear interpolation between 
E o (0) = -1.16944 J and E (l) = 0. 
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FIG. 5. Ground state energy estimated by high tempera- 
ture series as a function of doping at J/t = 2.0. The reverse 
curvature for 8 0.45 is due to extrapolating the high tem- 
perature free energy through T c > into the low temperature 
phase separated state. The reverse curvature gives a nega- 
tive compressibility, indicating the uniform phase is unstable 
towards phase separation. 

The ground state chemical potential in this parameter 
range is the constant slope of Eq(5) with the value \i/t = 
— 1.16944 J/t. The chemical potential hits the bottom 
of the tight binding band at J/t = 3.4367 and as J/t 
is further reduced the gain in kinetic energy eventually 
limits the phase separated state to J/t ^1.2. In Fig. 
6 we compare Eq to F(T) in the limit J/t — » oo with 
S = 0.5. 
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FIG. 6. Comparison of the free energy F(T) as a function of 
temperature to the known ground state energy Eq for 8 = 0.5, 
J/t — > oo. The crossing of these two curves we interpret as a 
phase transition at T c w 0.25J into a phase separated state. 
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Comparing the extrapolated F(T) to Eq we see they 
tend to cross at T ss 0.25 J. Since F(T) must be less than 
Eq and a monotonic function of T this crossing cannot 
occur. We interpret the tendency to cross as a phase 
transition to phase separation with T c « 0.25 J. 

Calculations for the 2D t-J model currently give a 
wide range of minimum J/t values for the presence of 
phase separation. Minimum J/t values reported in the 
literature are 0,|, 0.5-0.6 §,|,§, and our result of 
1.2 The latter results are in qualitative agreement 
in that there is a minimum J/t > for phase separa- 
tion. The reasons for these differences are not clear at 
present. However, while statistical errors are well under 
control, systematic errors in ground state energy calcula- 
tions due to small cluster sizes are much more difficult to 
control. Calculations investigating phase separation in 
the 2D Hubbard model (?J find e(5) equal to a constant 
for a range of dopings near half filling for the 17 = 
tight binding model. This spurious indication of phase 
separation is due to finite size effects and is reduced for 
larger clusters. Resolving the different reported results 
for phase sepraration will probably require significantly 
larger cluster sizes. 

In conclusion, by using an analysis of the high tem- 
perature series for the free energy per hole f(S) at differ- 
ent values of J/t we find that phase sepration in the t-J 
model is limited to J/t > 1.2. In addition, we find by in- 
direct arguments that T c ~ 0.25 J for the first order phase 
transition into the phase separated state. Combining this 
with the demonstration that phase separation can only 
occur at T — for the 2D Hubbard model on a square 
lattice |ll| supports the conjecture that the 2D Hubbard 
model does not phase separate for any positive U . Our 
results suggest phase separation in the 2D t-J model is 
a classical phase transition similar to a lattice gas with 
an attractive interaction [ fl2| and that phase separation 
is not important for physical choices of the t-J model 
parameters. 
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